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We present a derivation of the effective action for the relative phase of driven, aperture-coupled 
reservoirs of weakly-interacting condensed bosons from a (3+l)-D microscopic model with local 
U(l) gauge symmetry. We show that inclusion of local chemical potential and driving velocity 
fields as a gauge field allows derivation of the hydrodynamic equations of motion for the driven 
macroscopic phase differences across simple aperture arrays. For a single aperture, the current-phase 
equation for driven flow contains sinusoidal, linear, and current-bias contributions. We compute the 
renormalization group (RG) beta function of the periodic potential in the effective action for small 
tunneling amplitudes and use this to analyze the temperature dependence of the low-energy current- 
phase relation, with application to the transition from linear to sinusoidal current-phase behavior 
observed in experiments by Hoskinson et alP for liquid 4 He driven through nanoaperture arrays. 
Extension of the microscopic theory to a two-aperture array shows that interference between the 
microscopic tunneling contributions for individual apertures leads to an effective coupling between 
apertures which amplifies the Josephson oscillations in the array. The resulting multi-aperture 
current-phase equations are found to be equivalent to a set of equations for coupled pendula, with 
microscopically derived couplings. 
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I. INTRODUCTION 



A. Background 



The observation of coherent Josephson oscillations be- 
tween driven reservoirs of liquid 4 He separated by an 
array of nanometer-sized apertures^ has stimulated re- 
search into thermal/quantum fluctuations of the macro- 
scopic phase in aoerture array geometries^!. These stud- 
ies provide thermodynamic justification for the transi- 
tion, observed in Ref. pQ, between two different current- 
phase relationships as a function of temperature below 
the lambda point: the linear regime (occurring at low 
temperatures, T\ — T > 5mK ) in which the current de- 
pends linearly on the phase, and the "weak-link" Joseph- 
son flow regime in which the current has sinusoidal phase 
dependence (occurring for T\ — T < 0.8mK). The linear 
current-phase relationship at low temperatures is thought 
to be due to independent phase slips occurring at in- 
dividual apertures in the array in response to external 
driving by the hydrodynamic resonator. As the tempera- 
ture is increased toward the lambda point, the coherence 
length increases and the phase differences across indi- 
vidual apertures appear to become synchronized. It has 
been proposed that this results in coherent dissipative 
events, i.e. "phase-slip avalanches," giving way eventu- 
ally to coherent Josephson flow and a characteristic sinu- 
soidal current-phase relatiorP^. 

No microscopic quantum mechanical explanation cur- 
rently exists for this phenomenon. In order to justify 
the observed synchronization of the phase differences, 
Pekker, et al. postulate an effective long-range interac- 
tion between local phase gradients^ while Chui, et al. ex- 
ploit the analogy between a Josephson junction array and 



classical coupled pendula to explore thermal phase fluc- 
tuations in an aperture array^-l In this work, we derive an 
effective theory and equations of motion for the phase dif- 
ference across a single aperture and a simple two-aperture 
array, starting from a local U(l) gauge theory coupled to 
bosonic matter. The gauge field is necessitated by the 
presence of an external driving velocity which induces a 
"vector potential" v(r,t) and concomitant local chemical 
potential (f>(r,t), analogous to the electromagnetic gauge 
field in the theory of Bardeen, Cooper and Schrieffer 
(BCS) for superconductivity. 

Since our aim is a microscopic derivation of the equa- 
tions of motion for the macroscopic phase differences 



across aperture arrays (Section III I and an examination 



of the low-energy properties of the resulting current- 
phase relation across an aperture (Section [v]), we em- 
ploy here functional integral techniques rather than 
well known mean field or hydrodynamic techniques for 
bosonic systems (e.g., a gauged Gross-Pitaevskii (GP) 
cquatiorP for weakly-interacting Bose gases or a gauged 
two-fluid modeP for 4 He near the lambda point). The 
functional integral approach allows the equations of mo- 
tion to be derived from the microscopic Lagrangian, as 
was demonstrated for the analogous case of supercon- 
ducting systems by Ambegaokar et al. in RefP^ (hence- 
forth referred to as AES). Both the stationary phase anal- 
yses and the perturbative renormalization group proce- 
dure in this work are most convenient to carry out using 
this formalism. 

The fundamental variable of our effective theory is 
a gauge-invariant phase difference across an aperture: 
Aj(t) = A9(t) + m J dr ■ v(r,t), where the integral is 
taken on a short line segment through the aperture. It 
contains contributions from the background phase tex- 
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ture V0(r, r) (the irrotational superfluid velocity) and 
the external driving velocity (the gauge field). We show 
that the action governing the gauge-invariant phase dif- 
ferences for simple aperture arrays provides a qualitative 
explanation for the experimental observations of 4 He flow 
through nanoapcrture arrays over a range of tempera- 
tures below TxfM 

As noted above, a related microscopic derivation for 
a superconducting system appears in AES, in which an 
effective theory is derived for the dynamics of a super- 
conducting tunnel junction in terms of the macroscopic 
phase difference across the junction or the magnetic flux 
threading a superconducting quantum interference de- 
vice. Like AES, we shall be concerned here only with the 
dynamics of the low-energy degree of freedom in the sys- 
tem, namely, the macroscopic phase difference between 
junction-coupled Bose gases. In the present work, we 
focus on incorporating an externally imposed driving ve- 
locity into a gauge-invariant description of coupled reser- 
voirs of weakly-interacting bosons, on determining the 
current-phase relation for this system in different param- 
eter regimes dependent on the energy scale, and on using 
the results of this analysis to interpret the experimental 
observations of RefPl We shall not undertake the further 
analysis of real-time current correlations, dissipation due 
to quasiparticles, or the effects of noise in the junction 
that was also made in AES. Explicit comparison between 
our results for driven, weakly-interacting bosons with the 
results of AES for superconducting systems will be given 
where relevant in the subsequent sections. 



B. Summary of results 

The microscopic analysis presented in this work shows 
that the main features of the transition from linear-to- 
sinusoidal Josephson flow as a function of temperature 
are apparent already in the one and two-aperture cases. 
Starting from a local U(l) gauge-invariant Lagrangian, 
we derive the effective action for one and two-aperture ar- 
rays. We first show that a perturbative expansion of the 
gauge theory can be used to derive the quantum hydrody- 
namical equations of motion for the driven superfluid. In 
particular, we show that the Josephson-Anderson equa- 
tion for phase evolutiorP in gauge-invariant form, the 
circulation (superfluid fluxoid) quantization in the pres- 
ence of a driving velocity field, and the London equation 
leading to the Hess-Fairbank effect 9 can all be derived 
from the stationary-phase approximation to the effective 
action. A Legendre transformation of the Euclidean ef- 
fective action is then used to derive the current-phase 
relations for one and two-aperture arrays. We show that 
for a single aperture, the current-phase relation is consis- 
tent with a potential composed of sinusoidal, linear, and 
quadratic terms, while for a two-aperture array we find 
that interference between the microscopic tunneling con- 
tributions for individual apertures leads to a coupling of 
the current-phase equations of the two-aperture system. 



For the single aperture case, we then employ a weak- 
coupling renormalization group calculation to demon- 
strate the existence of temperature intervals in which the 
current-phase relation has predominantly linear or pre- 
dominantly Josephson (sinusoidal) behavior. The criti- 
cal temperatures separating these regions of linear and 
sinusoidal behavior are determined by relating the ratio 
of two coefficients in the rescaled effective action, each 
of which we calculate microscopically to one-loop order 
in perturbation theory, to the finite-temperature healing 
length. Application of the theory to the experiment in 
RefpJ on driven He flow through arrays of nanometer- 
sized apertures provides a rationalization for the transi- 
tion between linear and sinusoidal current-phase relation- 
ships that was observed as the temperature was increased 
toward the lambda point. 



C. Outline 

In Section |n] we discuss the local U(l) gauge invariant 
Euclidean action used in the coherent-state functional in- 
tegral and transform this action into a bilinear form in 
the real density field which can be analyzed using pertur- 
bation theory. In Section [TTT| integration over the density 
field is performed and the resulting perturbation series 
for the full inverse Green's function is used to determine 
the effective action for the phase difference across a sin- 
gle aperture. We show that the stationary phase approx- 
imation to the perturbed action allows gauge invariant 
forms of several superfluid hydrodynamical equations to 
be derived, e.g. the Josephson-Anderson equation for 
phase-difference evolution, the London equation for the 
gauge-invariant velocity, and circulation (superfluid flux- 
oid) quantization. The central result of this paper is the 
derivation of current-phase relationships for the single 



aperture and two-aperture array in Section [TV] We an- 
alyze the temperature dependence of the current-phase 
relation for a single aperture in the limit of small tunnel- 
ing amplitude by computing the RG beta function of the 
coupling constant Ej of the periodic potential and use 
this to analyze the experimental measurements of driven 
4 He flow through arrays of nanometer-sized apertures. 
We summarize in Section |VT| and discuss potential direc- 
tions for future research. 



II. THE MODEL 

We seek an effective theory for condensed, driven, 
weakly-interacting bosons separated by an array of one or 
two apertures in terms of local phase differences across 
the apertures. While our model includes only a local 
two-body potential, we will show that the main features 
of recent experimental results for liquid 4 He flow through 
nanoaperture arrays^ are nevertheless already explained 
by the current analysis. Our starting point is the Hamil- 
tonian in Eq. ([!]) for the weakly-interacting Bose gas that 
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is minimally coupled to a local chemical potential field s-wave scattering length), and D = V + imv(r,T) is the 

4>(r, t) and a vector field v(r, r) which will be interpreted covariant derivative. H ex t[v, <p] is a classical energy anal- 

as an external driving velocity. The Hamiltonian (with- ogous to electromagnetic field energy in superconductors 

out a tunneling term) is and depends only on external fields. 

Hl^ iZ>l = - / d 3 rDib 1< (r)Dtb(r) + — / d 3 r iii^ 2 (r)'ii) 2 (r) 

2m J 2 J We then construct the coherent state path-integral La- 

/grangian, Eq. (J2j) , corresponding to this Hamiltonian and 
d 3 r cj)(r)$ (r)ip(r) + H cxt [v, <f>] (1) additionally incorporate a single aperture tunneling term 

T r y that couples points r and r' on different sides of the 
The weak interaction is given by the usual delta function aperture. In the bosonic coherent state path integral, the 
two-body potential, with strength Vq (proportional to the Lagrangian is given by (H = kg = 1): 



Lhlj,ip*, A,V,d>] = I d 3 ri/j*(r,T)(d T + mcf>(r,T) - n)ip(r,T) + — [ d 3 rDip* (r, r) Dip (r, r) 

2m J 

+ / d 3 rd 3 r'i{;*(r,T)T r , r ,iP(r',T) + y J d 3 rA* (r, r)A(r, r) 
-;3 



d 6 r [A(r, r)if>* (r, r)ip(r, r) - A* (r, r)ip* (r, r)ip(r, r)] 



2 

"' L ~ f d 3 r(V x «(r, r)) 2 + ^ / d 3 r(i9 T w(r, r) - V0(r, r)) 2 (2) 



Here, A(r, r) and A*(r, r) are Hubbard-Stratonovich 
fields introduced to decouple the quartic interaction in 
the weakly interacting Bose gas, L has dimension of 
length, and r is the imaginary time. In the grand canon- 

- f* 3 drL 

ical partition function, Z(fi, j3) — J e Jo , the func- 
tional integration is over the fields ip, ip* , A, and A* 
and also the gauge field (with the measure defined in the 
discretized expression for the coherent state path inte- 
graP^). The last two terms are derived from H ext [v,(j)] 
in Eq.Q and are analogous to the electromagnetic field 
energy in superconductors; the vorticity (circulation en- 
ergy density) corresponding to the magnetic field energy 
density and an "electric" energy density analogous to the 
electric field energy density. The fields v(r, r) and 0(r, t) 
are analogues of the magnetic vector potential and local 
voltage of electrodynamics. These will be shown to sat- 
isfy stationary phase equations (Section III I and we do 
not analyze fluctuations of the gauge field configurations. 

If the tunneling matrix is multiplied by a U(l) parallel 
transporter via: 



im J dr-v(r,r) 



(3) 



this Lagrangian is clearly invariant under ip(r,r) — > 
ip(r, T)e lA ( r ' T ) (where A(r, r) is real) as long as v(r, t) — > 
w(r,r) - iVA(r,T) and (f>(r,r) <j>(r, r) - 4<9 r A(r,r). 
Put another way, we are analyzing a local U(l) gauge 
theory for the superfluid where (4>(r, t), v(r, r)) is the 



u(l) gauge field. The gauge transformation of the 0- 
component is due to working in imaginary-time (i.e. the 
base-space for the U(l) principal bundle is a Euclidean 
manifold). In the analysis to follow, it will lead to e.g. 
an imaginary Josephson-Anderson equation, which must 
be Wick rotated to obtain the real-time equation. The 
mean-field equations of the gauged weakly-interacting 
Bose gas are the stationary phase equations of this bare 
action: r ) = gives a gauged Gross-Pitaevskii 

equation^ for ip, while SA ^ r t) = =*> A(r,r) = 
-|V(r,r)| 2 and = A*(r,r) = \iP{r,r)\\ 

Note that the Hubbard-Stratonovich fields are not com- 
plex conjugates! This is a peculiarity of the bosonic 
Hubbard-Stratonovich transformation. Since the action 
does not depend on space-time derivatives of A(r, t) or 
A*(r, r), they may be taken as real constants at mean- 
field level. In the following, we choose A(r, r) = —A and 
A*(r, r) = A with A a real constant. 

To isolate a local phase field, a polar decomp osition 
can be made on ip and ip* , e.g. ip — » yj p(r, T)e ie ( r ' T \ 
This transformation does not change the measure in the 
functional integral for the partition function. Physically 
it means we are considering a restricted ensemble, i.e. we 
consider only a single condensed mode in the path inte- 
gral. This is our only explicit use of Bose-Einstein con- 
densation of the weakly-interacting Bose gas in this work. 
The action corresponding to the resulting Lagrangian can 
be brought into bilinear form: 



r0 



s = 



dr'dr 



d 3 r'd 3 r 



V^TjG-^r, r; r', r') ^JWy) + fA 2 S(r - r r )S(r - r') + — (V x v g (r, r)) 2 5(r - r')5{r - r') 
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+ n ^ L (idrv(r, t) - V0(r, r)) 2 5{r - r')S(r - r') 



(4) 



where v g (r,r) :— v(r,r) + ^V6*(r, r) is the gauge- 
invariant velocity. The operator G -1 (shown in Eqs. ^ 
- ([7])) is the object of principal computational interest in 
subsequent sections. Note that besides the field strength 
contributions and constant offset proportional to A 2 , the 
complete action can be written as a bilinear form. Us- 



ing the gauged GP equation^ and intepreting p( 
local condensate density field, it can be shown that the 
mean-field hydrodynamic effect of the external driving 
velocity is a depletion of condensate currenlP^K^ 



r,T) as a 



<9 T p(r, r) 



V.(-p(r,r)Vfl(r,r)) 
m 



V • (p(r, r)v(r, r)) 

In this article, we take the point of view that <p{r, t) and 
v[r, t) comprise the gauge- field in the fluid resulting from 
externally applied driving fields; in particular, the gauge- 
field is not internally generated by fluctuations. In the 
nanoaperture array experiment of Hoskinson et alP, the 
oscillations of the hydrodynamic resonator couple to both 
the condensate atoms and the depletion, like a piston. 
Thus, the gauge field can be viewed as the externally 
applied, nonconservative part of the total velocity of the 
fluid. The response of the phase field to the gauge field is 
apparent in the Euler and Josephson- Anderson equations 
that we derive below. An evolving velocity field induces a 
local chemical potential texture (via the Euler equation) 
which in turn induces an evolving phase field (via the 
Josephson- Anderson equation) . 



A. Self-consistent equation for A 

As mentioned in Section [TTJ since there are no space 
or time derivatives of the Hubbard-Stratonovich fields 
A(r, r), A*(r, r) in the action, we can take them to be 
constant. From the GP equation, |A| is equal to the 
density of condensed bosons and we take it to have the 
same value on both sides of the junction for simplicity. 
To compute the mean- field value of A from Eq. Q, we 
require that it extremizes the action: ^£ = 0. This 
mean-field is only present in G" 1 and in an additional 
term of quadratic order, and we use the Matsubara fre- 
quency and momentum representation of the free inverse 
Green's function to find the extremurrA^l Evaluating the 
resulting Matsubara surrf21 yields a self-consistent equa- 
tion for the mean-field A that is analogous to the BCS 
gap equation: 



A 



f 

V 



d 3 k 



(2tt) 3 e /9(J£-^+V A) _ 1 



(8) 



where V is the volume of the system. In deriving this 
equation we have assumed that the effect of the gauge 
field and tunnelling across the aperture contribute negli- 
gibly to the mean field value of A. 



B. G- , Josephson- Anderson equation 



III. PERTURBATION THEORY AND 
EFFECTIVE ACTION 

The operator G _1 (r, r; r', r') defining the bilinear form 
in the action Eq. Q is 

G" 1 = Gr 1 + G7 1 + G" 1 + G T l = Gr 1 + SG' 1 (6) 
where the individual components of <5G _1 are: 

Go" 1 = (d T - ^V 2 - M + V A)S(r' - r)S(r - r') 
G7 1 = [id T 6(r, t) + rruf>(r, r)] 5(r' - r)S(r - r') 

G~l = (\mv g (r,T) 2 )6(r' r)5(r r') 

G T l = T rr ,e^ r '^- 6i - r ^e im S r J ^(^ 6 (r - r') (7) 

Integration over the density field results in a term 
|tr logG -1 in the action. The trace is an integral over 
all internal positions or momenta and imaginary time 
t e [0, j3] arguments. Details of the perturbation ex- 
pansion for contributions to the action from each part of 
SG^ 1 are given in Appendi x A and the general techniques 
can be found in Reference d 6 * 14 * 15 !. 



In superconductors, the dynamical (a.c.) Josephson ef- 
fect is expressed by the Josephson- Anderson equation for 
phase evolution and is dependent on a voltage across the 
tunnel junction 16 . Although the weakly-interacting Bose 
gas is not charged, that does not preclude introduction 
of a 0-component of the gauge field, the fluctuating lo- 
cal chemical potential cj>(r, r), that appears in the action 
with the imaginary-time minimal coup ling. The first or- 
der contribution of <5G -1 = GJ 1 in Eq. ( Al I vanishes due 



to the periodic boundary conditions of (J{r, r) on [0, j3] 
and the requirement that the integral over both sides of 
the aperture (i.e. both reservoirs) is zero. However, a 
stationary phase equation with respect to (j)(r, r) can be 
derived from the action by finding the extremum of the 
"electric field" energy density term in Eq. Q, resulting 



id T v(r, t) — V0(r, r) 



(9) 



When transformed into real time, this becomes a classi- 
cal Euler equation relating the acceleration of the driv- 
ing velocity to a chemical potential difference across the 
aperture^. 
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The second order term in Eq. ( Al ) gives a nonvanishing 
contribution: 

-\J* dTG o (0, t- 0, rf (id T 6(0, r) + m0(O, t)\ (10) 

where the tilde signifies a move to momentum space. 
Go(0, r; 0, r) =: n is the (r-independent) number of k = 
bosons so it can be pulled out of the integral, result- 
ing in the coefficient — ^ (taking into account also the 
factor of 1 multiplying the perturbation series). The 
imaginary-time Josephson-Anderson equation can then 
be derived at this order from c ~ s , s . = together with 

the global phase and chemical potential configurations 
for each reservoir, 9{r,r) = r/l (t), 0(r,r) = (/>r/l(t), 
yielding: 

idrO^r) = -m^ L (r) 

td T 9 R (r) = -rafe(r) (11) 

Subtracting these equations gives the usual form of the 
Josephson-Anderson equation for evolution of the phase 
difference A9(t) across a junction: 

id T A9{r) = -mA</»(r) (12) 

We can now use the mean- field Euler equation, Eq. ([9]), 
derived from the bare theory to write the Josephson- 
Anderson equation in gauge-invariant form. Noting that 
4>l{t) = 4>r(t) + i f r ^ dr ■ d T v{r,r), we define the gauge 
invariant phase- difference by 

A 7 (r) = 9 l {t) - 6 r {t) +m f L dr- v(r, r) (13) 



and an analog of the electric field by 



c 



dr ■ (— V</>(r, r) + id T v(r, t)) 



(14) 



A rearrangement of Eq.(12)-(13) then yields the desired 



gauge invariant form of the Josephson-Anderson equation 

as 



id T Aj(r) = — mC(r). 



(15) 



To assess the contribution to the effective action S e s, 



Eq. (10 1 should be expressed in terms of Aj(t). ^>l(t) 
can be eliminated from the action using the Euler equa- 
tion and <Pr(t) can be eliminated by a Gaussian integra- 
tion (or vice versa, see All. The result is a capacitive 



term in the effective action: 



S, 



c 



drE c (d T Aj(T] 



(16) 



where the microscopic expression for Eq (from the Gaus- 
sian integral) is " |- . Since in our analysis the Eu- 
ler equation is considered a hard constraint, the elec- 
tric energy density vanishes and Sc is the only contri- 
bution from Gj 1 . Here and in other parts of the single 
aperture calculations, we neglect cross terms of the form 
tr log [GqSG^GoSG^ 1 ■ ■ •] which are, however, necessary 
for generating interactions between apertures in the mul- 
tiaperture case. 



C. G v 1 , circulation quantization, Hess-Fairbank 
equation 



Using Eq. (All to expand the contribution of G v ^ from 
Eq. ([7]) to first order results in a term quadratic in the 
gauge invariant velocity field V g (r, r) (i.e., a massive term 
for Vg). In superconductors, the physical consequence of 
a massive vector field is the Meissner effect, a repulsion 
of magnetic fields from the interior of the superconduc- 
tor up to a certain penetration depth which is dependent 
on the superfluid density^. The analogous effect for 4 He 
is the Hess-Fairbank effect, in which the superfluid mass 
density does not respond to rotation of the container due 
to an energy barrier to vorticity entering the superfluid 9 . 
The massive term for v g in the action suppresses fluctua- 
tions of the magnitude of the macroscopic phase gradient 
^ V# from that of the driving velocity v(r, r) in the bulk 
of the system. A stationary phase analysis of the action 



with respect to v g (r, t) at this order (see Appendix A 2 ) 



yields: 



J -v g {r,T)-L 2 V\(r,T) = 0, 



(17) 



which is a London equation describing the decay of the 
gauge- invariant velocity v g (r, r) in the interior of the 

bosonic system, with penetration depth A = \J 

Because we have included an external velocity field, it 
is useful to explore the consequences of this on circula- 
tion quantization. In the absence of driving (v(r, r) =0), 
the circulation integral is quantized in values of the cir- 
culation flux $0 = f?) due to the single- valuedness of 
the phase: f dr ■ v g (r, t) = ±§ V0(r, r) = I e Z. It 
seems clear that some form of the quantization should 
carry over to the driven case. To this end, we will 
integrate the London equation over a properly chosen 
contour. By analogy with Ampere's law, we interpret 
V x ui(t,t) as a current j (r, r) . We take a line integral 
of the London equation (171 around the torus on a path 



C which goes all the way around the torus, except for a 
missing segment C across the aperture (Figure [TJ 

dr - U(r,T) + ^(v(r,T) + ^S70(r,T)) \ =0 (18) 



Since 9 is only defined mod 2tt, one must have / dr ■ 

Jc 

V6»(r,r) = 2vrZ - A9(t), I <E Z where A9 is the line 
integral of the phase over C (local phase difference across 
the aperture). The line integral of the driving velocity 
gives the external circulation $ 

dr-v(r,r)+ / v(r,r) ■ dr = / to ■ dS = $ (19) 
c Jc J A 

Combining these two equations gives the generalized cir- 
culation quantization condition: 

1 4r. j{ r,r)\ = ^- e (20) 



2tt 



() 



FIG. 1. The contours used to derive superfluid fiuxoid quan- 
tization. C" is threaded through the aperture in the septum. 



with £ G Z an d A~/(t) the gauge-invariant phase differ- 
ence in Eq.(13). $0 is the circulation quantum, $0 = — • 



At distances into the bulk superfluid greater than the 
penetration depth, j(r, r) = so that the above equa- 
tion reduces to an equation for the quantization of the 
circulation due to the driving velocity. The general 
form of Eq. ( 20 ) expresses quantization of the superfluid 
"fluxoid"^ which contains contributions from vorticity 



due to the driving current in addition to the superfluid 
circulation. 

It remains to determine the contribution of v g to the 
effective action for the gauge-invariant phase difference. 
Inclusion of the vorticity (circulation energy density) in 
the bare action Eq. Q and requiring that the London 
equation hold results in the cancellation of the first order 
contribution of by the circulation energy density (see 

It should be mentioned 



Appendix A 2 for derivation) 
that in deriving this cancellation, we ignore a topolog- 
ical surface term J av v g A u). In fact, had we included 
in the Higgs action source terms for v g (r, r) and w(r, r), 
parametrized the vortex current by an appropriate gauge 
field, and integrated out v g (r, r) , the effective theory for 
the phase texture and vorte x gauge field would be a BF 
topological field theorjlM2!l. j n this work, we do not con- 
sider explicitly the dynamics of vortices (but see discus- 
sion in Section 



VII 



The second order contribution of gives a nonva- 
nishing contribution to the effective action for A7 and 
simplifies to: 



16 



d 3 kd 3 k' nkUk' 



d ^tf { ' 9< - 9 ' T ^( fc - k ' -1, r)v g {£, T)v g (k' - k 



(21) 



This expression is a convolution in momentum variables, 
but can be approximated as local in momentum because 
rifc is exponentially suppressed for k ^ at low temper- 
atures. Because the quadratic term in v g gives rise to 
a linear term in A7 (Eq. (A5)) the term quartic in v g 



D. 



Gj, 1 , periodic potential 



results in a quadratic term for Aj(t): 



S q = -E q (£- 



2tt 1 



(22) 



where Eq = m 2 n 2 \\v\\ 2 L 4 $1/16. 

We note here that in the present analysis of driven 
bosonic flow through an aperture, G" 1 is strictly second 
order in the gauge- invariant velocity v g , while the corre- 
sponding perturbative contribution for superconducting 
current flow thr ough a Josephson junction also contains a 
ullSl The consequence is that our second 
is quartic in v a . This difference 
is a result of the polar decomposition of the bosonic fields 
made here into real components, in contrast to the su- 
perconducting case in which one must work with Nambu 
spinors. Expanding the square in Eq. (22 ) shows that the 



term linear in 
order expansion in G~ 



effective action has both quadratic and linear dependence 
on A7. The latter will result in a quantized constant term 
(a quantized current-bias) in the current-phase equation 
while the former will give a term proportional to A7 (see 
Section ITV]. 



We require that the tunneling matrix T rr > — when 
r and r' are on the same side of the aperture and, for 
simplicity, T rr > = T = const, when r and r' are on oppo- 
site sides of the aperture. In the perturbation expansion, 
we must integrate over all possible positions which give 
nonzero tunneling matrix elements. The resulting term 
in the effective action is: 



S.j = Tn 



f 

/ drcosA7(r) 
Jo 



(23) 



where n is the zero-momentum occupation. If the per- 
turbation expansion is continued and the imaginary-time 
integrations are approximated by a single one, higher har- 
monics of the cos(A7(r)) interaction result; we will not 
include these in our analysis. These interactions can be 
shown to be of less relevance than the leading interaction 
(decreasing faster as the high-energy cutoff is lowered) by 
background field RG methods^!. However, if one keeps 
imaginary-time arguments distinct (i.e., preserves time 
non-locality) in the second order contribution, the sec- 
ond order term may be included as a dis sipat ive contri- 
bution to the effective action (see Section A3). AES use 
an analogous term of this order to model the effect of 
quasiparticle-macroscopic phase difference scattering on 
the current in the Josephson junction. 
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IV. EFFECTIVE ACTION AND 
CURRENT-PHASE RELATIONS 



to be S c s = S c + S Q + Sj. Explicitly: 



A. Effective action 

The effective action for the gauge-invariant phase dif- 
ference is determined from Eq. (16 1, Eq. (22), Eq. (23) 



S bS [&7(t);1,P\= cLtE c 9 T A 7 (r) - E Q (t 



A 7 (r), 
2tt ' 



+ Ej cos A 7 (t) 



(24) 



where the microscopic expressions for the coefficients 
have been derived above: Eq — n 2 V 2 /8 ,Eq — 
m 2 n 2 \\v\\ 2 L 4 /16 ,Ej — Tn. This effective action de- 
scribes a particle on a ring with a potential that is a sum 
of a parabolic and cosine terms, i.e., 

V[Aj] = -E Q (£- A 7 (t)/2tt) 2 + Ej cos A 7 (r). (25) 

(see Figure [2]). In the partition function involving S e g, 
the sum over £ € Z counts the winding number of the 
macroscopic phase. Changes in I correspond to phase 
slips across the aperture. The behavior of the effective 
potential for £ — and a range of relative values of the 
parameters Ej,Eq is shown in Figure [2j Note that for a 
given £ there is an infinite number of local minima of the 
potential. The generalized circulation quantum condi- 
tion, Eq. ( 20 ) , may be used to further write the effective 
action solely in terms of the circulation The Hamil- 
tonian corresponding to this action is formally similar to 
that used to desc ribe rf SQUIDs and superconducting 
flux qubits^ 23 * 24 ! and has been used previously to ana- 
lyze coherent quantum phase slips^. 

The quadratic contribution of A7 in the effective po- 
tential differentiates this potential from the sinusoidal- 
plus-linear or "washboard" form of effective potential 
found for a current-biased Josephson junction^. The 
effective action derived here for driven bosonic flow 
through an aperture differs from that derived by AES 
for superconducting flow through a Josephson junction 
in two respects. First, for the driven bosonic flow, the 
gauge field contribution G^ 1 to the effective action at 
second order in perturbation theory is quartic in v g , re- 
sulting in a term quadratic in A7 and hence a parabolic 
contribution to the potential. In contrast, the contribu- 
tion from the superconducting superfluid velocity to the 
effective action for a superconducting tunnel junction is 
linear in the phase difference variable (see Eq. (31) in 
RefP) and second order terms arise only from the addi- 
tional inductive energy. Second, we have neglected the 
second order tunneling perturbation which is nonlocal in 
time: inclusion of this would, by analogy with the analy- 
sis of AES, give rise to dissipation in the aperture array. 

In Section [V] below we will analyze the temperature- 
dependence of Ej. Because the temperature-dependence 



will enter through the ratio of Eq to Eq in Eq. ( 24 ) , we 
now show that the latter ratio can be written in terms 
of the ratio of two characteristic lengths of the system. 
According to the analysis above: 



Eq = 27r 2 \\v\\ 2 L 4 
E c V 2 

^ 2 \\v\\ 2 n 2 X 2 L 
2V 2 



(26) 



where \l is the penetration depth of the gauge invariant 
velocity v g 



(Eq. 



17). The condensate density, n/V, is 
related to the healmg length at nonzero T in the Popov 
theory by n/V = l/87ra£ (T) 2 , where a is the s-wave scat- 
tering lengtrPS We can then express the ratio by 



E c {2aV2) 2 V2£(T) 



(27) 



We will use Eq. (271 to analyze the current-phase relation 
of a single aperture in Section [IVB 1| below. 



B. Generalized current-phase relations 

The current-phase relation is the central equation of 
weak link Josephson phenomena. Not only does it reflect 
the macroscopic quantum nature of the flow through an 
aperture, but it can also yield indirect information on 
the microscopic dynamics of the constituent particles at 
the aperture. The most extensive studies of such flow 
for interacting bosons in atomic systems (as opposed to 
superconductors) have been made for 4 He, where exper- 
iments with driven flow through nanoaperture arrays re- 
veal the existence of two different current-phase relations 
in different temperature regimes below the A point. We 
will use the microscopically derived effective action for 
the single aperture, Eq. (24) to construct the current- 



phase relation for the weakly interacting bosonic system 
and use it to analyze the 4 He flow experiment, bearing 
in mind that liquid 4 He is a strongly interacting system 
so the analysis remains qualitative. Since observation of 
Josephson effects under external driving has not yet been 
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FIG. 2. The potential of the effective Euclidean action, 
V[A7] = —Eq(£ — A7/271-) 2 + _EjcosA7 versus A7 shown 
with I = for E.j = 1 and three values of the ratio Eq/Ej: 
Eq/Ej = 10 (blue), 1 (green), 0.1 (red). Evaluating the 
potential at i 7^ breaks the A7 — > — A7 symmetry. 



observed for weakly interacting Bose condensed gases, al- 
though both Josephson coupling 27 and persistent flow 28 
have been observed in different geometries (double well 
and toroidal traps, respectively), we also expect that our 
analysis will be applicable to driven condensed weakly- 
interacting Bose gases separated by aperture arrays. 



Current-phase equation for single aperture 



The current-phase equation resulting from the effec- 
tive action S e s is obtained as the stationary phase equa- 
tion SSes/SA^yir) — 0. This equation is derived in con- 
venient form by first defining the "density difference" 
field An(r) that is canonically conjugate to A7(r), by 
Legendre transformation of the Lagrangian in the path- 
integral. Specifically, the kinetic term of Eq. (24) is 
changed via: 



9xA 7 (t) 



V[An(r)]e 



- dr ^An(rf+i ^ dr An(r)3 T A 7 (r) 



(28) 



Performing a stationary phase analysis with respect to 
A7(t) on the resulting Legendre transformed Eq. (24) 



then yields the general imaginary time current-phase 
equation 



E, 



id T An(T) - Ej sin(A7(r)) + — ^-A7(t 

27T^ 



E, 



Q' 



:((29) 



The current-phase relation in terms of the real-time cur- 
rent I(t) = dAn/dt is obtained by a Wick rotation of 
Eq. ( 29 ) . The term linear in A7 confirms that this 



current-phase relation constitutes an analog for weakly 
interacting condensed bosons of the generalized Joseph- 
son equation for an rf SQUID. The current-bias part 
of the current-phase relation is constant and quantized, 
proportional to I 6 Z. We emphasize that Eq. (29) con- 



tains all terms necessary to describe a linear-to-sinusoidal 
current-phase transition. 

The different forms of the current-phase relation in 
different physical regimes correspond to specific values 
of the parameters Eq and Ej. For Ej — 0, the 



current-phase relationship of Eq. ( 29 1 is linear and cor- 



responds to the small amplitude oscillations of a pendu- 
luir l 29 | 30 [ However, the effect of an ^-dependent current- 
bias persists. For Eq = 0, this equation reduces to the 
imaginary-time version of the Josephson equation, with 
critical number current equal to Ej — Tn (and mass 



current given by mEj). 

To determine (in imaginary time) the classical equa- 
tion for A7, one can require the exponent of Eq. (28) 
to be stationary with respect to variations in An. This 
results in the relation idrA^ = ^j-An, analogous to 
mq = p in classical mechanics. Substituting this relation 
into Eq.(29), it is then evident that for Eq <C Eq the 
quantized current-bias and the coefficient of the linear 
term are negligible; it is in this regime that purely sinu- 
soidal oscillations should be observed. In this limit, one 
recovers the imaginary time version of the classical (fixed 
length) pendulum equation with amplitude Ej /2Eq, i.e., 



9 2 A 7 



Ej 
2E C 



sin(A7), 



(30) 



which constitutes a well-known classical analogue of the 
Josephson effeclP^. 



More generally, the current-phase relation, Eq. ( 29 1 , 



interpolates between two regimes of purely linear and 
sinusoidal current-phase equations at Ej = and Eq = 
0, respectively (plotted in real time in Figure [3]). These 
two limiting current-phase behaviors were observed for 
different temperature intervals in the 4 He nanoaperture 
array experiments of RefPl 

If we set = A ap in Eq. (27), with A op the di- 
ameter of a single aperture (see Section [V] for justi- 



9 




FIG. 3. Plots of real time (i£ -> -f t in Eq. (291 ) 
current-phase relations for mixed linear-sinusoidal contribu- 
tions over [— 7T, it] for £ = and a range of values for the 
parameter ratio Ej/Eq at fixed temperature. These are to 
be compared with the experimentally observed current-phase 
relations shown in Figure 3 of Reference^. In Section \V\ we 
derive the temperature-dependence of the current-phase rela- 
tionship. 



fication), we may relate the ratio Eq/Eq to the ra- 
tio A ap /2£(T) of the characteristic aperture size to the 
temperature-dependent healing length. As the temper- 
ature is decreased, the healing length becomes smaller 
than the aperture size and the ratio Eq / Eq grows quar- 
tically. Thus if Ej is considered fixed, the linear term 
in the current-phase equation Eq. ( |29[ ) becomes dom- 
inant for low T. In contrast, at higher temperatures, 
e.g. large enough that the healing length is larger than 
the aperture size, the sinusoidal term would become 
dominant. Whether thermal fluctuations of the gauge- 
invariant phase difference wash out the sinusoidal part 
of the current-phase relation as T\ is approached from 
below depends on the size of Ej, the scaling of which is 
derived in terms of E c and Eq in Section [V] 



This qualitative analysis shows that as the temperature 
is increased towards T\, there can be a transition from 
a linear current-phase relation at low temperatures to a 
sinusoidal current-phase relation at higher temperatures 
(but still below T A ). 



2. Current-phase relation for two- aperture array 



Within the framework of this theoretical analysis, 
adding an additional aperture is straightforward and re- 
sults in a substantially richer set of current-phase phe- 
nomena. We analyze here just the two-aperture case, 
leaving the extension to arrays with large numbers of 
apertures for future investigation. We may assume the 
cross-sectional areas of the two apertures are identical. 
There are now two tunneling matrices T^, and T^,,; we 

require that T^, is nonzero only when r and r' are on 
opposite sides of aperture 1 and both are in a small vicin- 

(2) 

ity of the aperture (similarly for J, ) . In addition to the 
sum of single aperture effective actions for the gauge in- 
variant phase differences A^ x \t) and A^ 2 \t), which 
have been derived in Section |IIID| there is now also a 
tunneling cross-term that appears at second order in the 
perturbation theory. This tunneling cross-term generates 
an effective aperture interaction that may be expressed 
in terms of the microscopic phase differences across the 
individual apertures. In particular, with the tunneling 
amplitudes assumed to be the same, this term adds an 
interaction to the effective action for two apertures of the 
form 



5 lnt = -£}cos(A 7 W(t))cos(A 7 ( 2 )(t)). 



(31) 



For small phase-differences, expansion of this equation 
implies that the homogeneous part of degree 2 renormal- 
izes the quadratic parts of the uncoupled contributions to 
the action and introduces a coupling A7W (t)A^ 2 \t) , 
while the homogeneous part of degree 4 introduces a cou- 
pling A7^ 1 ^(r) 2 A7^ 2 - ) (r) 2 as well as quartic local poten- 
tials for the phase differences. Neglecting these higher 
order terms, the interaction results in coupled modi- 
fied Josephson equations which describe classical coupled 
pendula. 



We can use the two-aperture coupling term Eq. (31 1 to 



rationalize the experimentally observed transition from 
a linear to sinusoidal current-phase relation in a multi- 
aperture array. Because the coefficient of the interac- 
tion just derived is the square of Ej we know that if Ej 
is large compared to Eq, the current-phase relation for 
each individual aperture is approximately sinusoidal and 
that the energy cost for having an inter-aperture phase 
difference of tt is 2Ej. This means that for Ej ^ 0, it is 
favorable for the difference of the phase-differences to be 
mod 2-7T. Hence the amplitude of the oscillation coming 
from the independent terms is doubled. This is consis- 
tent with both the experimental observations of phase 
difference synchronization as the current-phase relation 
becomes sinusoidal, i.e. Josephson-like, as well as with 
the observed linear scaling of the Josephson oscillation 
amplitude with number of apertures^. 
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V. RENORMALIZATION GROUP ANALYSIS 
FOR SMALL E,, 

To make contact with experiment and to justify the 
qualitative argument presented in Section |IVB 1[ it is 
desirable to understand how the current-phase relation- 
ship of the effective theory, Eq. ( 29 ) , and in particular 
the critical current Ej, depends on temperature. This 
can be done by employing RG methods in the small Ej 
regime and analyzing the corresponding beta functiorP^. 
The sign of this function determines how the coupling 
constant Ej behaves (i.e. decreases or increases) at low 
energies/long length scales. 

Full details of the RG calculations are included in Ap- 
pendix B. Here we summarize only the key features of 
this calculation and the results that are relevant to un- 
derstanding the temperature dependence of the current - 



2ir(2irE' c + 1) 



phase relation presented in Section IV B above. We note 
that in order for the system to be described by the phase- 
difference only, we must implicitly assume a high-energy 
cutoff A, beyond which energy scale the effective theory 
is invalid. At the energy scale determined by b = 4, 
with A a lower energy scale (i.e., b e [l,oo)), Ej(b) is 
the critical current of the current-phase relation and its 
magnitude relative to Eq will determine the Josephson 
character of the current-phase relation. 

Since we are concerned here with the scaling of Ej (b) , 
we neglect the scaling of Ec and Eq. If Ej{b) decreases 
(increases) as we consider low energy scales, we infer that 
the low-energy current phase relation Eq. p9| does not 
contain (does contain) a sinusoidal term. The resulting 
beta function is then given by 



dEj 
1 db 



16=1 



1 



2ttA 



Att 2 E c A 2 - E c 



Ej. (32) 



Integrating this differential equation by separating vari- 
ables and transforming to dimensionless parameters (us- 
ing the naive scaling dimension of each) E'q = 
E' c = EgA, yields the following scaling field for Ej: 



E'c 




FIG. 4. Renormalization group diagram for the Josephson 
(sinusoidal) contribution to the generalized current-phase re- 
lation, Eq. ( |29| |. The solid line separating regions I and II 
is the marginal line E'q = 2n (2nE' c + 1) and the solid line 
separating regions II and III is the singular line E'q — A-r 2 E' c . 
A plus (minus) sign in a given region indicates the direction 
of the Ej RG flow, corresponding to the cosine potential be- 
ing relevant and increasing with increasing 6 (irrelevant and 
so decreasing with increasing b). A relevant periodic poten- 
tial results in a Josephson term in the current-phase relation. 
Each temperature T defines a ray in this positive quadrant 
(see text). 



Eq/Ec considered at the beginning of Section IV (see 
Eq. i27|)). In order to evaluate this ratio as a function of 



the renormalization scaling b, we must choose a value for 
the high-energy cutoff, A. In the low temperature helium 
nanoaperture array experiments of RefP, the largest en- 
ergy scale is the kinetic energy of the driving velocity. 
We therefore employ an energy cutoff value A = 2 ^jL . 
Returning to Eq. (27 1, we see that the condition for the 



singularity will then occur at a temperature T\ such that 



Ej(b) cxEjb (1+s \ 
1 



( X i 



(33) 



We have confirmed the validity of this scaling field with a 
background RG calculation^. When the exponent 1 + g 
is negative, Ej(b) will be irrelevant and disappear at low 
energies, while when the exponent is positive Ej(b) is rel- 
evant and grows at low energies. Figure [4] shows the re- 
sulting RG flow diagram for Ej in the positive (E'q, E'q) 
quadrant. 

There are two important features in this RG diagram 
for Ej. First, the singular line defined by E'q = ^t 2 E' c 
(where the denominator of g goes to zero) and second, the 
marginal line at E'q = 2tt (2nE' c + 1) (where 1 + 9 = 0). 
We can analyze the singular line in terms of the ratio 



E'c \^{Ti) 



4tt^ 



(34) 



In the following analysis of the experiment in RefP} 
we set Xl = A ap , with A ap the diameter of an aperture 
in the array. This is consistent with analysis of the first 
critical angular velocity for appearance of vortex lines 
in rotating annular reservoirs of liquid 4 He, where the 
annular width appears in the expression for critical an- 
gular velocity in an identical form to the London pen- 
etration depth for first critical magnetic field in type- 
II superconductors^^!] From the known temperature- 
dependence of the healing length £(T) for He iPH, and 
using the experimental aperture diameter of 40 nm, this 
scaling singularity is found to occur at a critical temper- 
ature Ti w T\ — 20mK. Conversely, for a general tem- 
perature T, the right-hand side of Eq. (34 1 shows that 
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each value of T defines a ray in the coupling-constant 
space (E' c ,E'q), and specifying the exact point in cou- 
pling constant space requires knowledge of either E' c or 

We now consider the nature of the current-phase re- 
lation for temperatures below and above T\ . Figure [4] 
summarizes the structure of the low-energy current-phase 
relation in the (E'q,E' c ) plane for the regime of small 
Josephson coupling Ej. At temperatures below T-y, 
E'q/E' c > Air 2 and we are either in region I or region II 
of Figure [4] In the former case we might expect a mixed 
sinusoidal/linear flow, while in the latter case we expect 
only a linear current-phase relation. In region I, E'q is 
always nonzero so some linear flow is always present. Al- 
though the experiment in RefP does not address this par- 
ticular temperature regime, we can use a number of ar- 
guments to predict the expected balance between linear 
and sinusoidal contributions as a function of temperature 
within this regime. At low temperature when the number 
of condensed bosons is large, or whenever the tunneling 
amplitude is very large or very small, the Ej cos A7 part 
of the action can be treated using the Villain approxima- 
tion^ which would renormalize Sq and lead to a purely 
linear current /phase equation. In support of this argu- 
ment is the fact that for Eg > 2ir and constant, the 
value of E' c is lower in region I than in region III. A 
low value of E' c implies a high energetic cost for density 
difference fluctuations (see Eq. (28). Since the density 



difference is canonically conjugate to the gauge-invariant 
phase difference, we expect that a low variance in the 
value of the former quantity allows for a high variance in 
A7 and hence for the Josephson flow contribution to the 
current-phase relation to be washed out. 

At temperatures above T\, the periodic potential is 
relevant (region III). The current-phase relationship, 
Eq. ( 29 1 will always have a nonvanishing contribution 



from sinusoidal flow in this regime (while the system re- 
mains below T\ although it may be mixed with linear 
flow. For E'q small, nearly pure Josephson oscillations 
should be observed. 

To support the validity of this analysis of the small 
Ej current-phase relation, we place two results from the 
experiments of RefP that exhibit different current-phase 
behaviors into the context of the RG diagram, Figure [4] 
For example, at T\ — T — 27mK a linear current-phase 
relation is observed. Employing the experimental for- 
mula for the healing length^ and an aperture width 
A Q p = 40nm, yields the ray E'q/E' c « 94 for this tem- 
perature. Since the experimental current-phase relation 
has linear character at this temperature, we expect that 
this point lies in region II below the Eq = 2tt(2ttEq + 1) 
line. The second point we analyze is T\ — T = 0.8 mK. 
Here the experiment shows nearly pure Josephson oscil- 
lations and experimental estimates for healing length and 
aperture width yield the ray E'q/E' c s» 8.0 x 10~ 3 . Con- 
sequently this higher temperature point lies in region III, 
far below the E'q = E' c line and in a region where E'q is 
negligible. 



We emphasize that pure sinusoidal Josephson oscilla- 
tions (without the modified dynamics due to parabolic 
potential) should be found in region III of Fig. ([I]) only. 
This is a regime of considerable interest for applications 
of Josephson phenomena in liquid 4 He to metrologj^fil 
and for development of cir culation analogues of super- 
conducting flux qubitiP^^. The experimental challenge 
in accessing this regime lies in the fabrication of small 
enough nanoaperture arrays in order for the E'q = 4ir 2 E' c 
line to be reached deep in the condensed phase and not 
near the critical point. 



A. Multiple apertures 



The present analysis is made for a single aperture. Ob- 
serving a Josephson current for a bosonic superfluid in a 
single driven nanoaperture is known to be a challenging 
task, due to the small amplitude of oscillation compared 
to the amplitude of oscillations of the driving device. Our 
analysis shows that if the healing length of an interact- 
ing Bose gas can be made over twice the characteristic 
aperture size, nearly pure Josephson oscillations would be 
observable. Unfortunately, for driven liquid 4 He in aper- 
ture arrays of X ap ~ 40nm, the system for which all such 
experiments have been performed to date, this regime is 
nearly precluded by the lambda transition. In supcrflu- 
ids with larger zero-temperature coherence lengths (e.g., 
the paired fermion superfluids, including 3 He and many 
type-II superconductors) the Josephson effect is conse- 
quently more robust with a single aperture. 

For bosonic superfluids such as liquid 4 He and trapped 
dilute Bose gases, it is of interest to consider what 
changes to the present analysis are required by having 
multiple apertures. If tunneling amplitudes at each aper- 
ture are the same and each aperture has the same size 
and shape, even the particulars of the weak Ej coupling 
RG calculation should carry over. The most important 
change in going from one aperture to multiple apertures 
is the presence of the phase-difference interaction and 
the independent tunneling terms as mentioned in Sec- 
If the phase-difference interaction favors a 
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tion 

uniform value, the classical configurations will be phase- 
locked, independent tunneling terms will add up and the 
overall tunneling amplitude will be scaled by M, with 
M the number of apertures in the array. Consequently, 
the amplitude of the Josephson oscillation is multiplied 
by M and it is easier to observe. It should be noted 
that the presence of multiple apertures introduces new, 
higher-order operators in the effective action. In general, 
their anomalous scaling dimensions (and hence their op- 
erator relevance) are different from that of the cos(A7) 
potential. 
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VI. CONCLUSION 

In this work we have derived and analyzed an effec- 
tive theory of gauge-invariant phase differences across 
simple aperture arrays starting from a local U(l) gauge 
theory. The stationary-phase approximation to the lo- 
cal U(l) gauge theory at first and second order expan- 
sion of the one-loop contribution to the action was shown 
to reproduce many well-known equations of motion, e.g. 
the Josephson-Anderson equation, the Euler equation, 
the London equation, the equation of superfluid fluxoid 
quantization, and the d.c. Josephson equation. We have 
shown that the general current-phase relationship is con- 
sistent with the phase dynamics in a potential formally 
analogous to that of a rf SQUID, consisting of quadratic, 
linear and sinusoidal terms whose relative strength is de- 
termined by the magnitudes of the charging and Joseph- 
son couplings, Eq and Ej, respectively. The effective 
action leading to this current-phase relationship differs 
from that derived by AES in the context of supercon- 
ductive tunneling^ due to the explicit presence of the 
parabolic potential in the action, as well as to the locality 
in time assumed in our analysis. Analysis of dissipation 
in the aperture array deriving from the second order time 
nonlocal contribution of G^ 1 will be addressed in future 
work. 

The effect of the sinusoidal term in the current-phase 
relation was further analyzed using finite temperature 
renormalization group methods. We have shown that 
the sinusoidal part of the current-phase relationship is ex- 
pected to become significant in two different regimes, but 
that it is most important when the coherence length £(T) 
is larger than the characteristic size of the aperture, X ap . 
By exploiting the relationship between Eq/Ec and the 
ratio of the aperture size to the temperature-dependent 
healing length, we were able to examine the scaling of 
Ej with respect to this ratio. This analysis identified 
regions II and III, separated by a singular line in the RG 
diagram, that are respectively consistent with the linear 
and sinusoidal current-phase relations that were observed 
experimentally in RefP. Using the relevant experimental 
values of healing length and aperture dimensions, we have 
shown that the singular line separating these regions, 
Eg = 4tt 2 E' c , occurs about 20 mK below the lambda 
transition. The qualitative agreement of this value with 
the experimentally observed transition at ~ 5 mK below 
T\ in Ref. lj provides strong evidence for the validity of 
this effective theory. In addition, generalization of the 
effective action derived in this theory from one to two 
apertures shows that phase-difference coupling between 
multiple apertures leads to phase-difference synchroniza- 
tion and to a doubled amplitude of Josephson oscillation 
in the array. Our analysis indicates that for M parallel 
apertures in an array, we may expect the amplitude of 
Josephson oscillations to behave as O(M). 

In this paper we have considered neither the dynamics 
of phase slips and the vortices by which they are carried, 
nor their role in the transition from linear to sinusoidal 



current-phase relationship (see Section IIICl. However, 
we note that inclusion of the nonlocal interaction be- 
tween uj(r) and ui{r') is expected to lead to the hydrody- 
namic equations first presented in reference |39j . In the 
nanoaperture array, the low-temperature linear current- 
phase characteristic is thought to be due to independent 
nucleation and subsequent slippage of vortices at individ- 
ual apertures^. These events dissipate the kinetic energy 
of the hydrodynamic resonator slowly, as opposed to large 
scale coherent phase slips occurring at higher tempera- 
tures. In this regime, the diameters of vortex cores are 
nearly as large as the apertures themselves. This suggests 
a physical picture of vortex proliferation at the nanoaper- 
ture array leading to coherent oscillations. Such a picture 
is consistent with our requirements that i) Ej be relevant 
in order to observe Josephson oscillations, and ii) Eq be 
small so that the gauge-invariant phase difference is not 
pinned to an integer multiple of 2n. In the core of vor- 
tices pinned at the array, off-diagonal long range order is 
destroyed and A7 is allowed to fluctuate away from 27r£. 
The dynamics of the vortices may be studied by deriving 
their effective theory using boson- vortex duality 40 . Such 
a study would be useful both to confirm in the dual pic- 
ture the features of the phase diagram derived here, and 
to investigate the properties of a vortex condensate in an 
aperture array for which the bosonic field operator used 
here no longer describes particles above the vacuum. 

In utilizing the current approach to interpret exper- 
iments on liquid He, we have neglected the strongly- 
interacting nature of superfluid helium, i.e., we do not 
consider a realistic two-body potential. Realistic stud- 
ies for Josephson effects in liquid helium driven through 
nanoscale aperture arrays may be undertaken with path- 
integral Monte Carlo methods^. To our knowledge 
there has so far been no observation of Josephson os- 
cillations between driven reservoirs of weakly-interacting 
condensed bosons separated by nanoaperture arrays, nor 
indeed of any Josephson effects under driving flow con- 
ditions for weakly interacting Bose condensate systems. 
However, the Josephson effect has bee n obs erved for 
weakly coupled Bose-Einstein condensate d 27 * 42 *, and per- 
sistent currents have been observed in toroidally trapped 
condensates^. Taken together with the recently demon- 
strated ability to make arbitrary potentials in such ge- 
ometries^, the rapid progress in experimental study and 
manipulation of rotating BECs in toroidal traps holds 
out the prospect of future realization of Josephson phe- 
nomena in confined atomic BECs. 



Appendix A: Perturbative expansion of G 



For convenience and clarity, the perturbative expan- 
sion of G _1 is included in this appendix. We use the fol- 
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lowing perturbative series to analyze the action Eq. Q: 



tr log [Go 1 + SO' 1 ] = tr log [G " 1 ]+^ 



\fe+i 



-tr [(Go^G" 1 ) 
(Al) 



The first term in this series is a constant which can- 
cels due to the normalization of the partition function. 
The second term, a series in powers of <5G _1 , gives im- 
portant contributions to the effective action. The free 
Green's function of the action is found by inverting the 
Go 1 operato^ 



G (k, r; k', t') = (2n) 3 5{k - k') exp[-£ fc (r - r')](6(r - r')(l + n k ) + 0(r' - r)n k ) 



(A2) 



where £t 



*1 

2»i 



(j, + VqA and n k 



1 



is the 



Bose- Einstein distribution. When r = r', the time- 
ordered correlation function is the normal ordered cor- 
relation function and so Go(k, r; k , t) = rik5(k — k). 

In evaluating the integrations over internal momenta, 
we frequently use the fact that Go(k, r; k', t') oc S(k-k'). 
Treating the perturbation series exactly results in nonlo- 
cal contributions to the action. We assume when needed 
that the imaginary time arguments of the higher order 
terms are the same, by appealing to the fact that the 
free Green's function is exponentially suppressed as dis- 
tance in imaginary time increases. Momentum integrals 
over the free Green's function are restricted to k — be- 
cause we are considering the low-energy dynamics of the 
condensed mode. 



G: contribution 



The tr log expansion with respect to th is perturba- 
tion is outlined in the text (Section III E> I . To convert 
to a functional of A7 in the effective action, 
is split (in position space) into left and right 



Eq 
Eq 



parts as 



2 _Yl 
4 



/(f dr Ud T 9 L + m<t> L ) +(L R). Take 

Vl = Vr = V for simplicity. Using the Euler equation, 
4>l is eliminated from the action. We then perform the 
Gaussian integral over <f> R to arrive at an effective term 
involving Aj(t) only. The Gaussian integral is: 



T>[4> R ] exp 



dTV 2 n 2 m 2 (j) 2 R 



oc exp 



Vn 2 mi f p ( f VL 
— - — J dr [drOh + d T 9 R + m J dr ■ d T v(r, r) ) (j) R 



2t/2 rP 



n 2 V 



J dr (^d T ^L + d T 6 R + m J dr ■ d T v(r, r)^ 



(A3) 



The (^^-independent part of the contribution is added 
into the exponent and the square expanded, yielding Sc 
after simplification. Elimination of the 0-component of 
the gauge field from the action is reasonable because it 
is not a dynamical field. 



2. G„ contribution 

Analyzing the contribution of the vorticity energy den- 
sity and the first order term in G,^ 1 to the perturbation 
expansion in Eq. ( Al I shows that the gauge-invariant ve- 
locity satisfies a London equation. We will use this equa- 
te 



bution to the action is 



tion to show that the first-order contribution of G„, 
the effective action for A7 is canceled by the vorticity 
energy density term. Specifically, the first-order contri- 



1 

2 tr 

nm 



G (k,T;k,T) ■ —v g (q,T)v g (-q,T) 
dr j d 3 r v g (r, r) 2 , 



where n is the number of condensed bosons. We have re- 
stricted the sum over momenta in the free Green's func- 
tion to k — because contributions from k 7^ are ex- 
ponentially suppressed at low temper atur es. 



To derive the London equation Eq. 

at first order in the expansion in G 
the Euler equation, the identity V 



we set 



ss_ 



= 



and make use of 
x (V x v g (r,T)) = 
t) , and the physical requirement 



V(V-Vg(r,T))~V 2 v g (r, 
that v(r,r) be divergence-free. Owczarek has exploited 
a similar "Higgs"-type argument to rationalize the ex- 
pulsion of circulation by a superfluid, noting that if the 
source of v g [r, t) is a roton, the penetration depth is 
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roughly the same as experimentally-observed vortex core 
diameters!^ 



In Section [1110 
tribution of G" 1 



we stated that the first order con- 
is canceled by the circulation energy 
density. This can be seen as follows: consider integrat- 
ing the first order contribution along an integral curve T 
of the superfluid velocity. Assuming a toroidal geometry 
as in the derivation of superfluid fluxoid quantization in 
Section |III C[ we can approximate the first order contri- 
bution of v g by 



nm 



d 3 r v g (r, tY 



nmL 



dsv g (r,T) 2 , (A4) 



density, which can be rewritten 
mL 2 



d r v g (r,T)j(r,T) 



(A6) 



The vector identity a ■ (V x b) — b ■ (V x a) — V • (a x 
b) has been used in deriving this formula. Substituting 
into Eq.(A6) the London equation in the form j(r,r) = 
^jj vqjr , t), one obtains the perturbation contribution in 
Eq. ( A5 ) but multiplied by a factor of — 1 . In using the 



vector identity above, we have neglected a topological 
contribution to the effective action. This is discussed in 
Section UTTCl 



contribution 



where the latter integral is with respect to arc-length 
and L 2 is the area factor multiplying the vorticity energy 
density in the microscopic Lagrangian. This integral can 
be converted to a line integral by identifying the tangent 
vector to T with the driving velocity at each point. This 
is justifiable because i) F is an integral curve of the su- 
perfluid velocity and it is physically reasonable to assume 
that for low T, V#(r, t) is parallel to v(r, r) at each point 
in spacetime, and also since ii) ||V0(r, r)|| <C \\v(r, t)||. 
For a constant magnitude driving velocity, the first order 
contribution becomes: 



nmL 2 f 2 77J7j||v||L 2 f 

— ^— J dsv g (r,T) = j^dr-Vg{r,T) 



4 Z7T 



This contribution is canceled by the circulation energy 



In this calculation, as in previous ones, we specialize 
to the left/right reservoir phase configuration 9(r,r) — 
9r/l( t )- In the multiaperture case, these become lo- 
cal left /right macroscopic phases in the vicinity of each 
aperture. Employing the convention that the left-to- 
right gauge-invariant phase difference is defined to be 
— A7(t)), the first order contribution is 



1 



,tr 



rtr 



Go(r L ,T;r R ,T)T rRrL e 



zA7(t 



G a {r' R ,T-r' L: T)T r> ,e^^ 



The corresponding first-order contribution Sj to the 
effective action results from using the fact that 
J d 3 rd 3 r'G {r, r; r' , r) =G (k-k' = 0, r). 

In our analysis of the effective theory for the gauge- 
invariant phase difference, nonlocal imaginary time terms 
in the perturbation expansion have been neglected. Here 
we derive one of these nonlocal terms arising from the 
second order contribution of ; diagrams correspond- 
ing to this contribution are shown in Fig. 

The resulting contribution is: 



exp 



\ ftr [G (r L , r; r R , ^T^e^^G^y '; r' R , T)T r , RrL e-* A ^ + (R o L)] 



+ tr [G (r L , r; r' Ll r')T r , r , e iA ^G (r' R , r'; r R , r)T rRrL e- 1 ^ + (R o L)] 



(A7) 



Transforming to momentum space and taking the tun- 
neling amplitude to be a constant, T, yields: 



T 
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:)(,r-r') 



(A8) 



^ 2 J drdr' J d 3 kd 3 k'e^ 

(1 + n k )n k , cos(A7(r)) cos(A7(t')) 

Our expression for the nonlocal contribution for this 
driven bosonic flow differs from that of AES because we 
do not have a particle-hole symmetry. 



Appendix B: RG for periodic potential 



We start by expressing the effective action Eq. ( 24 ) in 
terms of Matsubara frequencies (we use the S c g label for 
both the imaginary-time and Matsubara representations 
of the action): 



S eS [A T ,l,P] 



n=0 



E c < 



Eq 
4tt 2 



A7(w„)A7(w_ 
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• l — r R 



Go 

Grr, 



r' T ^ 













FIG. 5. Diagrammatic representatio n of the two second 



order contributions from G T l in Eq.(A7) which result in 



an imaginary time nonlocal contribution to the effective ac- 
tion. The convention — A7 is used for left-to-right hopping. 
We omit the delta function vertices as they are omitted in 
Eq.dA7l. 



dr Ej cos A7(r) 



f^A 7 (0).(Bl) 



Choosing a high-energy cutoff A, A 7 (r) can then be 
split into low-frequency (slow, s) and high-frequency 



(fast, /) terms, A7 s (r) 



*T/(r) = /A <H<A &e-~ 



t A 7 (cj), respectively, where 
b is the renormalization scaling. The effective action is 
then split into slow (s), fast (/) and combination (U) 
components: 

S cS [A7; l,0\ = S s [Aj] + S f [Ay] + Su[A 1s (t) + A 7/ (r)] 
with Su[f(T)] = / drSjcos[/(r)]. (B2) 



We note that the slow part gets an additional contribu- 
tion from the zero mode in Eq. (Bl). Assuming that 



T <C A, so that the Matsubara sums become integrals, it 
is then possible to integrate over the fast components by 
making use of a small Ej approximation^ 



-Sl„ 



]=e- s -[A7d (1 _ Su [A 7s ,A 7/ ] + ...) / 



e -S 5 [A 78 ] e -(S c/ [A 7s ,A 7/ ]> / 



to obtain an effective low energy action SWwenergy [A 7a ]. 
Here the / subscript denotes an expectation value using 
Sf as the action. 

An explicit evaluation of (SuiAjs, Ajf])f results in a 
6-dependent multiplicative renormalization of Ej, which 
we call Ej(b). This integration over fast modes is given 
explicitly by: 



(Sc[A 7s) A 7/ ]) f=Ej I Z>[A 7 (w)]e 
><|w|<A 



C.C. [B3) 



Carrying out the Gaussian integration (and neglecting quency still lies within the positive interval [0, 00) results 
the divergent contributions) gives in Eq. ( 32 ) in the main text and we see that the peri- 

odic potential is multiplicatively renormalized. We have 
verified our result for the /3-function using a background- 
field RG analysis according to the procedure outlined in 
Ref.DE]. 



(Su) f =Ej I dre 



' cos(A 7s (t)) 



(B4) 

Rescaling r —¥ ? to ensure that the Matsubara fre- 
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